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Abstract. Simultaneous observations of the accretion-powered millisecond pulsar IGR J00291+5934 by 
Inter-national Gamma-Ray Astrophysics Laboratory and Rossi X-ray Timing Explorer during the 2004 December 
outburst are analysed. The average spectrum is well described by thermal Comptonization with an electron tem- 
perature of 50 keV and Thomson optical depth tt 1 in a slab geometry. The spectral shape is almost constant 
during the outburst. We detect a spin-up of the pulsar with v = 8.4 x 10~^^Hz s~^. The ISGRI data reveal the 
pulsation of X-rays at a period of 1.67 milliseconds up to ~ 150 keV. The pulsed fraction is shown to increase from 
6 per cent at 6 keV to 12-20 per cent at 100 keV. This is naturally explained by the action of the Doppler effect 
the exponentially cutoff Comptonization spectrum from the hot spot. The nearly sinusoidal pulses show soft lags 
with complex energy dependence, increasing up to 7 keV, then decreasing to 15 keV, and seemingly saturating at 
higher energies. 

Key words, accretion, accretion discs - binaries: close - pulsars: individual (IGR J0029H~5934) - stars: neutron 
- X-ray: binaries 



1. Introduction 

IGR J00291-I-5934 was discovered in the galactic disk 
by the Internati onal Gamma-Ray As trophysics Laboratory 
{INTEGRAL) l)Eckert et all l2004j) during a routine 
Galactic Plane Scan of the Gas A region on December 

2, 2004 during an outburst. Follow-up observations on 
December 3, 2004 with the Rossi X-ray Timing Explorer 
{RXTE) classified the source as an accreting X-ray mil- 
liseco nd pulsar (MSP), the s ixth known source of this 
class l)Markwardt et alJl2004a|) . IGR 100291+5934 has a 
spin period of 1.67 ms, making it the fastest known X- 
ray MSP to date, and an orbital period of 2.5 hours 
l)Markwardt et alJl2004b|l . The other known X-ray accret- 
ing millisecond pulsars' spin frequencies lie between 180 
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and 435 Hz (see review bv lWiinandjl20n5l) . Their orbital 
periods fall into two distinct ranges - either around 40 
minutes or 2-4.5 hours. 

A candidate optical counterpart was observed with di- 
minishing intensity from magnitude i? 17 — 21, consis- 
tent with the X-ray flux decay rate of IGRJ00291+5934 
during the latter days of the outburst iIFox fc Kulkarnil 
l2004t iBikmaev et al.H200l ISteeehs et al.ll2004jl . It is the 
second X-ray MSP (after SAX J1808.4-3658) for which 
the optical flux decay could be observed. Spectroscopy of 
the optical source le d to the detection of broad emission 
lines of HcII and Ha (jRoelofs et al.l200^lFilipDenko et alJ 
l2004b Radio observations also detected a source consis- 
tent with the posit i on of IGR J002914-5 934 l)Poolevll20o4 
iFender et aLll2004t iR^pen et al J 12004*). Th e source has 
been observed in quiescence with Chandra bv i.Tonker et alJ 
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and a first INTEGRAL spectrum was reported by 
IShaw et alJ ll200fj^ . 

The hard X-ray emission of IGR J00291+5934 was 
significantly detected by INTEGRAL during the entire 
fifteen-day observation from outburst to quiescence. In 
this paper we analyze the X-ray spectrum during the out- 
burst. Furthermore, we perform timing analysis on the 
INTEGRAL high energy data as well as on data from a 
RXTE follow-up observation, in order to study the char- 
acteristics of the pulse profile (shape, time lags, pulsed 
fraction) at energies from 2 to 150 keV. 



2. Observations and data 

2.1. INTEGRAL 

The dataset was obtained with INTEGRAL 
ijWinkler et al.l 120031 from December 2 to December 16 
2004, i.e. from satellite revolution 261 to 265 (see Table 
nj. These revolutions include a Target of Opportunity 
(ToO) observation and part of a Gas A/Tycho A02 
observation. 

We use data from the code d mask imager IBIS/ISGRI 
l)Ubertini et alJl2003l: iLebrun e t al. 200j) at energies from 
20 to ~ 200 keV (total exposure of 437.4 ks) and from the 
JEM-X monitor l)Lund et al.ll2003|) at energies 5 to ~ 20 
keV (total exposure time of 194.6 ks). For ISGRI, the data 
were extracted for all pointings with a source position off- 
set < 12°, and for JEM-X with an offset < 3.5°. The spec- 
trometer (SPI) was not used since the source was too weak 
for the SPI and its 2° angular resolution is inadequate in 
this case (see section IS^ . Data reduction was performed 
using the standard Offline Science Analysis (OSA) version 
4.2 distributed by the I NTEGRAL Science Data Genter 
llGourvoisier et al .120031) . The algorithms used for the spa- 
tial an d spectral analysis are described in lGoldwurm et all 

(Eool. 
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2.2. RXTE 

After the initial follow-up observations by RXTE start- 
ing December 3, 2004, a ToO observation was performed 
between December 7 and 21 (observation id. 90425), for 
which the data were made publicly available. For our tim- 
ing analysis, we also used these data, namely from the two 
non-imaging X-ray instr uments, the Proport ional Gounter 
Array (PGA; 2-60 keV) l|Jahoda et alJll9 96^ and the High 
Energy X-ray Timing E xperiment (HEXTE; 15-250 keV) 
jRothschild et a,lJh998ll . 

The PGA data were collected in the E_125us_64M_0_ls 
event mode, recording event arrival times with 125 ^s time 
resolution, and sorting events in 64 PHA channels. Default 
selection criteria were applied and the final net PGU2 ex- 
posure was 105.5 ks, for the December 7-10 data. We in- 
cluded in our analysis HEXTE ON-source data, using de- 
faul screening criteria for Gluster and 1. The screened 
dead-time corrected HEXTE Gluster and 1 exposures 
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Fig. 1. The 20-40 keV ISGRI mosaicked and deconvolved 
sky image of the ^ 437 ks observation. Image size is 
~ 6?25 X 3°, centered at IGR J00192-H5934 position. The 
pixel size is 5'. IGR J00291+5934 and V709 Gas were de- 
tected at a significance of ~ 88a and ~ IScr, respectively. 

were 35.7 ks and 36.4 ks, respectively, for the December 
7-10 data. 



3. Results 

3.1. IBIS/ISGRI image and light curves 

Fig. n shows a significance map in the 20-40 keV en- 
ergy range centered on IGR J00192-I-5934 located in the 
Gas A/Tycho region. Single pointings were deconvolved 
and analyzed separately, and then combined in mosaic 
images. Two sources are clearly detected at a signifi- 
cance level of 88.4cr for IGR J00291-I-5934 and 18.2a 
for the nearby intermediate polar V709 Gas. In the en- 
ergy band 40-80 keV, the significance level was 51.2a for 
IGR J00291-I-5934 and 6a for V709 Gas. At the higher en- 
ergies 80-200 keV the confidence level dropped to 17.1ct 
for IGR J00291-H5934, while V709 Gas was not detected at 
a statistically significant level neither in single exposures 
nor in the combined image. To obtain precise source lo- 
cations we simultaneously fitted the ISGRI point spread 
function to the two close sources. We obtained a posi- 
tion for IGR J00291+5934 at aj2ooo = 00^29'"02!92 and 
(5.12000 = 59°34'06'.'4. The position of V709 Gas is given 
by aj2ooo = 00'^28™55';29 and 5j2ooo = 59°16'14'.'0. The 
source pos ition offsets with respect to the optical cat alog 
positions ||Fox fc Kulkarnil2004t iDownes et alJll997j) are 
0'.2 for IGR J00291-h5934 and 1.'5 for V709 Gas. The er- 
rors are 0'.2 and 1.'5 for IGR J00291-I-5934 and V709 Gas, 
respectively. These are within the 90% confi dence level 
assuming the source location error given by iGros et alJ 
(2003). The derived angular distance between the two 
sources is ~ 18'. Due to the fact that INTEGRAL is able 
to image the sky at high angular resolution (12' for ISGRI 
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Table 1. Log of INTEGRAL observations analysed in this paper. Start and end times are in days of December 2004 
(from MJD 53340). Exposures are in seconds and count rates are in counts per second, for ISGRI in the 20-100 keV 
and for JEM-X in the 5-20 keV energy bands. 
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Fig. 2. INTEGRAL/ISGRI light curve in the 20-100 
keV energy band (averaged over 0.1 day intervals). The 
INTEGRAL data from the whole observation (Table [U 
have been converted to flux assuming a Comptonization 
model (see Table E)). The optical flux decay in R-mag is 
show n with diamonds , wher e the first two points arc taken 
from lFox fc Kulkarnil l)2004j) and lSteeehs et al.lil200l. re- 
specti vely, and the rest was taken from iBikmaev et all 
l|2004l) . The dashed hnes correspond to F oc e"*/^'^"* and 
Foce-*/2-2". 



and 3' for JEM-X), we were able to clearly distinguish 
and isolate the high-energy fluxes from the two sources 
separately. This allowed us to study the X-ray emission of 
IGR J00291+5934 during its entire outburst. 

The INTEGRAL 20-100 keV high energy light curve 
has been extracted from the images using all available 
pointings and is shown in Fig. [21 averaged over 0.1 day 
intervals. T he first point corre sponds to the detection of 
the source l|Eckert et al.ll2004j) . Using a Comptonization 
model (see section 15. 2 .in . we estimated the peak bolomet- 
ric X-ray flux to be 2.1 x 10~^ erg cm~^ s~^. For a distance 
of 5 kpc (see section E!T|l , this corresponds to a bolometric 
luminosity of 6.3 x 10'^® erg s~^, or 3.5% of the Eddington 
luminosity, LEdd) for a 1.4 Mq neutron star. This is a 
lower value then measured for XTE J1751-305 which ha s 
a luminosity of ~ 0.13LEdd ((Gierlinski fc Poutanenl2005l) . 



but similar to that of the other accreting X-ray MSPs. 



The outburst profile of IGR J00291-f5934 (Fig. EJ 
shows a deca y similar to those of four other millisecond 
pulsa rs (e.g., Gilfanov et al.lll998t IGierlinski fc PoutanenI 
l2005j) . Only XTE J1807-294 shows a pu rely exponen- 
tial de cay with a time scale of ^^120 days l|Falanga et alJ 
I 200,'j) . After the peak the flux declines exponentially, with 
a decay time-scale of 6.6 days (10 days in SAX J1808.4- 
3658, 7.2 days for XTE J1751-305), until it reaches a 
break, after which the flux drops suddenly with a de- 
cay time-scale of 2.2 days (1.3 days in SAX J1808.4- 
3658, - 0.6 days for XTE J1751-305). The higher energy 
(ISGRI) decay time scale is fully consistent with the JEM- 
X time scale at lower energy (3-20 keV). The decay of the 
optical flux, also shown in Fig. [3 is similar to that of the 
X-ray light curve, but somewhat smoother, with a decay 
time-scale of 4.4 days. Note that the optical source was 
visible up to 45 days after the discovery. The ISGRI light 
curve was also extracted in the following energy ranges: 
20-40 keV, 40-80 keV and 80-200 keV. The hardness ra- 
tio of these energy bands as a function of time indicates 
that no significant spectral variability was detected. 

3.2. Spectral analysis 

We verified that during the whole observation V709 Cas 
was observed up to 80 keV with a constant mean count 
rate of ~ 0.6 in the 20-80 keV energy band, and does not 
influence our spectral analysis for IGR J00291-I-5934. We 

J ierformed the spectral analysis using XSPEC version 11.3 
Arnaudl [l996l). combining the 20-200 keV ISGRI data 
with the simultaneous 5-20 keV JEM-X data. A constant 
factor was included in the fit to take into account the 
uncertainty in the cross-calibration of the instruments. A 
systematic error of 2% was applied to the JEM-X/ISGRI 
spectra which corresponds to the current uncertainty in 
the response matrix. All spectral uncertainties in the re- 
sults are given at a 90% confidence level for a single pa- 
rameter (Ax^ = 2.71). 



3.2.1. Total spectrum 

The joint JEM-X/ISGRI (5-200 keV) spectrum was 
first fitted with a simple model consisting of a 



4 



M. Falanga et al.: INTEGRAL observation of IGR J00291+5934 



photoelectrically-absorbed power-law (pl). Given that we 
were not able to constrain the iVn value (as the JEM-X 
bandpass starts above 5 keV) we fixed it to t he value found 
from Chandra observations at lower energies (jNowak et al.l 
l2nn4li . A simple PL model is found to be inadequate with 
X^/dof=105/36. The addition of a cut-off significantly im- 
proves the fit to x^/dof=47/35, with the bes t-fit photon 
index of ~ 1.8 (similar to that reported bv IShaw et al.l 
120011^1 and a cut-off energy of ~ 130 keV. However, this 
model does not describe the spectrum well below 10 keV, 
which is more complex. 

The emission from the accretion shock on a 
neutron star is expected to be produced by ther- 
mal Comp tonization of soft s e ed photons from the 
star (e.g. Iz^ldovich Shakural 1196^ lAlme fc Wilso^ 



Table 2. Spectral parameters of the thermal 
Comptonization COMPPS fit to the JEM-X/ISGRI data 



EomllCierlinski fc Pouta,nenll20oi lFa]a,nga et al.ll 

We model the shock emission by the COMPPS model^ 
l)Poutan cn & Svcnsson 1996), where the exact numerical 
solution of the Comptonization problem in different ge- 
ometries is obtained by considering successive scatter- 
ing orders. We approximate the accretion shock geometry 
by a plane-parallel slab at the neutron star surface. The 
main model parameters are the Thomson optical depth tt 
across the slab, the electron temperature Tq, and the soft 
seed photon temperature Tsood- The emitted spectrum de- 
pends also on the angle between the normal and the line 
of sight 9 which does not coincide with the inclination 
of the system because of light bending. The seed pho- 
tons are injected from the bottom of the slab. A fraction 
exp(— Tt/ cos 61) of these photons reaches the observer di- 
rectly, while the remaining part is scattered in the hot 
gas. Thus, the total spectrum contains unscattered black 
body photons and a hard Comptonized tail. The best fit 
with x^/dof=44/37 is then obtained for fcTgccd « 1.5 keV, 
fcTe « 50 keV, TT ~ 1.1, and cos6' « 0.6 (see TableE)). The 
model also allows us to determine the apparent area of the 
seed photons, which turns out to be ^sccd 21(Z?/5 kpc)^ 
km^ in the best fit. This corresponds to a hot spot radius 
of ^ 2.5 km in this outburst phase. The unabsorbed EFe 
spectrum and the residuals of the data to the model are 
shown in Fig. |31 

Soft emission (presumable from the accretion disk) 
with peak temperature of about 0.5 keV has been 
found in other accreting X-ray MSPs (se e e.g. 
iGierlihski fc PoutanenI boOSt iFalanga et alJ I 20051) using 
XMM-Newton data in the 0.5-10 keV energy band. As 
the JEM-X bandpass begins at 5 keV, it is impossible 
to search for this emission in these data. Neither the 6.4 
keV iron line nor Compton reflection were significantly 
detected. 



kTo (keV) 
fcT.ecd (keV) 
Tt 

^s^ccd (km^) 
cos 9 
xVdof 

T a 

^1-30C 



0.28 (f) 

491^ 
1.49j 



-1+0.16 
-0.32 
1 19+0.04 
^■^^-0.07 
12.6 
4.5 
+0.06 



(10^'' erg s" 



20.7 

0.60_o.o9 
44/37 
3.7 



° Assuming a distance of 5 kpc. 



Il973t ILvubarskii fc SunvaevI Il982) . Indeed, spectra from 
accretio n-powered MSP s are w ell described by this 
model llCierlihski et al.l I 2002t IPoutanen fc Gierlihskil 
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Fig. 3. The unfolded spectrum of IGR J00291+5934 fitted 
with an absorbed COMPPS model. The data points corre- 
spond to the JEM-X (5-20 keV) and ISGRI (20-200 keV) 
spectra, respectively. The total spectrum of the model is 
shown by a solid curve. The lower panel presents the resid- 
uals between the data and the model. 

3.2.2. Outburst spectrum 

We analysed the one day averaged JEM-X/ISGRI spectra 
for the observation during the outburst (see Table We 
plot the best fit results using the thermal Comptonization 
COMPPS model in Figure 0] The inclination angle, 6, and 
the seed photon temperature, Tgood, were fixed at the best 
found fit value for the entire observation (see Table . 
The luminosity Lboi was calculated for a distance of 5 kpc 
from the best fit model in the energy range 1-300 keV. 

The results show that the decay of the outburst is 
marked by a nearly constant plasma temperature, except 
for the last day of outburst where it shows a decrease along 
with an increase of the scattering optical depth. Neither 
variation, however, is very significant statistically. 

3.2.3. Total fluence 

From the broadband spectral model we estimated a to- 
tal fluence from December 2-15 in the 0.1-300 keV en- 
ergy band to be 1.37 x lO""^ erg cm^^. The fluence esti- 
mated with RXTE is higher by ~ 0.4 x 10~^ erg cm~^ 
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Fig. 4. Tlie outburst evolution of the best-fit spectral 
parameters of the COMPPS model. Each point, except the 
last, corresponds to a one day averaged spectrum; due to 
the lower flux level the last points are averaged over two 
days. 
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Owing to the 1° field of view of 
RXTE, the flux difference could be due to the intermedi- 
ate polar V709 Gas, located 20' from IGR J00291-h5934. 
The RXTE pointings were stable around q;j2ooo — 
00'"29'"09?23 and (5j2ooo = 59°33'58'.'l, i.e. at every point- 
ing where both sources were in the RXTE field of view. 
We analysed the spectrum of V709 Gas with the identical 
JEM-X and ISGRI data and exposure times reported in 
Table ^ The source was detected from 5 to 80 keV. In 
order to estimate the bolometric flux during the outburst 
period, we fitted the spectrum of V709 Gas up to 80 keV 
with a thermal bremsstrahlung model. The best fit with 
plasma temperature ~ 20 keV gives a x^/dof = 12/15. 
The flux in the energy band 0.1-100 keV was found to 
be 2.1 X 10-1° erg cm'^ 8"^ The flux of V709 Gas was 
constant throughout the IGR J00291-I-5934 outburst, so 
the total fluence contributed by V709 Gas from December 
2 to December 16 was ~ 0.3 x 10"'^ erg cm~^. The dif- 
ference between our measured total fluence and that of 
RXTE can therefore be attributed in large part to the 
flux of V709 Gas. The remainder is likely due to differ- 
ences in the models used to estimate the source flux and 
uncertainties in the cross-calibration of the instruments. 



3.3. Timing characteristics 
3.3.1. The pulsar ephemeris 



Searching the low-statistics ISGRI {E > 20 keV) 
frequency space for the 1.67 ms signal from 
IGR J00291+5934, involves making numerous trials, 
which reduce substantially the sensitivity to periodic 
signals. Therefore, we have used the high-statistics PGA 



data (2-60 keV), in order to obtain an accurate ephemeris 
for the millisecond pulsar. 

We first determined the orbital parameters by deriving 
the instantaneous pulse frequency during sufficiently short 
exposures (as compared to the orbital period) of typically 
200 s, in which the pulse profile could still be significantly 
recognized. We apphed the Zj'-statistic (|Buccheri et alJ 
^983) to pulse phase distributions for trial frequencies 
in a small window centered on the expected pulse fre- 
quency. The PGA data until 2004 December 11 (included) 
yielded sufficiently high signals for ^ 200 s integration 
times. The derived distribution of frequency and frequency 
uncertainty versus barycentered time (mid of integration 
interval) was subsequently subjected to an epoch folding 
method yielding the "best", optimized orbital period Porb- 
Folding with this best period resulted in an orbital phase 
versus frequency distribution, which is highly sinusoidal 
without any significant asymmetry. The amplitude result- 
ing from a sine/cosine fit gave the value of the projected 
semi- major axis of the neutron star orbit sini, and the 
minimum of the fit yielded the time of the ascending node 

('^asc ) ■ 

The derived orbital characteristics based on the PGA 
data from December 7-11 are ful ly consistent with those 
reported bv lGallowav et alJ ||200,'tI) based on December 3-6 
data. Therefore, in the further a nalysis we will use the set 
of orbital parameters derived bv lGallowav et all 1 200fJ ^ in 
order to facilitate direct comparisons (note Tasc = ~ 
0.25 • Porb = 53345.1619261 MJD TDB). 

At this point, we corrected the barycentered time tags 
of the PGA events for acceleration effects along the or- 
bit. For each PGA sub-observation we determined the 
time of arrival (TOA) by a pulse profile template corre- 
lation analysis similar to the technique normally used in 
radio -pulsar studies (see e.g. Tempo, ^ Flhv lo r fc Wcisbcrgj 
Il989|) . The template arrival times could be determined 
accurately only for PGA data collected between 7 and 
12 December 2004 (MJD 53346.202-53351.967), when 
IGR J00291+5934 was sufficiently strong. Phase folding 
the template arrival times through a timing model with 
a c onstant frequ ency ly — 598.89213064(1) Hz, as found 
bv I Gallowav et al.l l|2005|) . gave rise to large systematic 
deviations up to 0.15 in phase from the predicted arrival 
times. A timing model with 2 parameters, v and re- 
moved these deviations (all the residuals, TO As minus 
model, were < 0.025 in phase). The fitted timing model 
has: V = 598.89213060(1) Hz; z> = +8.4(6) x lO'^^ jj^, 
s^i at epoch MJD 53346 (TDB). Note that we measure 
a significant spin-up during this outburst. An error in the 
assumed source coordinates can give rise to timing er- 
rors which may introduce a s purious spin-up or spin-down 
l|Manchester fc Tavloil Il977(l . The average 0.092 arcsec 
source position error given bv iRupen et all l|2004|) would 
introduce a non-existant spin-up rate of ;> = 5.8 x 10"^^ 
Hz s~i during our observation. Since this is smaller than 
our measured frequency derivative, we conclude that the 



http:/ /pulsar. princeton.edu/tempo 
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source was spinning up during the outburst. An appar- 
ent spin-up on the order of our measured v would require 
a fairly large ~ 0.7 arcsec source position error during 
our observation. The ephemeris and orbital parameters 
are given in Table |31 



Table 3. Ephemeris of IGR J00291+5934 



Parameter Value 

Right Ascension (J2000)l' 00"29™3':0822 

Declination (J2000)t 59°34'18'.'99 

Epoch validity start/end (MJD) 53341 - 53353 

Frequency 598.89213060(1) 

Frequency derivative 4-8.4(6) x 10~^^ 

Epoch of the period (MJD;TDB) 53346 

Orbital period 8844.092 s 

Oajsini 64.993 (It-ms) 

Eccentricity 

Longitude of periastron 0° 
Time of ascending node (MJD;TDB) 53345.1619261 



Significant pulsed emission was detected up to 100 
keV, with some indication for a signal also between 100 
and 150 keV. The Z\ significance for deviations from 
uniformity was: 9.1cr, 7.3(7, 5.0cr and 2. Oct for the 20-35 
, 35-60, 60-100 and 100-150 keV energy bands, respec- 
tively. The sinusoidally shaped pulse profiles are shown in 
the panels k-n of Fig. The 1.67 ms timing signal was 
also detected in data from the JEM-X detector: 4ct non- 
uniformity significance for the 5-10 keV band (no event 
-selections except on energy were applied; see Fig. panel 



To complement and verify our ISGRI timing results we 
Hz also produced pulse phase distributions based on HEXTE 
Hz s~^(15~250 keV) data from the same time period as the PCA 
events above. The HEXTE profiles are shown in panels 
g-j of Fig. The Z\ significances were: 11.4ct, 8.3ct, 3.3ct 
and I.Ict (=no detection) for the energy bands, 20.3-35.2, 
35.2-60.1, 60.1-100.9 and 100.9-151.1 keV, respectively 



^ The position is based on radio observations with the 
VLA bv Rupen et al. f 20041 . 



3.3.2. INTEGRAL/RXTE pu\5e profiles and time lags 

We corrected the time stamps of the PCA events from the 
2004 December 7-10 (MJD 53346.185 - 53350.015) ob- 
servations to arrival times at the solar system barycenter 
using the JPL DE200 solar system ephemeris and the po- 
sition of IGR J00291-h5934 given in Table 01 taking into 
account the binary nature of the system. The barycen- 
tered times were phase-folded using the pulsar ephemeris 
given in Table |21 This yielded pulse phase distributions 
in 64 energy channels. Pulsed emission was detected up 
to ~ 60 keV. The non-uniformity significance derived by 
applying a Zj-statistic is even 5.7ct for the 34.7-58.1 keV 
band. The pulse profile shapes remain highly sinusoidal 
over the ~ 2 — 60 keV PCA range (see Fig. ^panels a-e). 

The ISGRI data from INTEGRAL revolutions 261 - 
263 (see Table ^ were used in our timing analysis aim- 
ing at the detection of the 1.67-millisecond timing sig- 
nal. INTEGRAL science windows (sew) showing erratic 
(ISGRI) count rate variations, indicative of effects due to 
Earth radiation belt passage or solar fiare activity, were 
excluded from the analysis. Only time stai nps of events 
with rise times between 7 and 90 channels ("Lebrun et al.' 
Hool, detected in non-noisy ISGRI pixels which have an 
illumination factor of more than 25% were passed for fur- 
ther analysis. The event times of these selected events are 
barycentered, and the subsequent pulse phase folding of 
the barycentered event times using h', v and the epoch of 
the period (see Table OJ yielded pulse phase distributions 
for different energy bands between 20 and 300 keV. 



Note that the pulse profiles detected by the PCA, 
JEM-X, HEXTE and ISGRI are fully consis tent in shape 
and absolute timing. iGallowav et al.l l|2005l) reported an 
energy-dependent time delay, with the 6-9 keV pulses ar- 
riving up to 85/iS earlier than those at lower energies (see 
Fig. 2 bottom panel of Galloway et al. 2005). Exploiting 
the ISGRI and HEXTE pulse phase distributions above 
20 keV, we can investigate this behavior up to ^ 100 
keV. We performed a similar analysis, first determining 
the time lags between profiles in different PCA energy 
bands, but now for the independent PCA dataset from 
2004 December 7-10. We confirm the reported variation 
in time lag between ~ 2 — 20 keV as a function of energy 
(Iflallowav et al. 20051) and extend the PCA coverage up 
to - 35 keV. 



Correlating the three significant ISGRI and two sig- 
nificant HEXTE pulse profiles between 20 keV and 100 
keV with the low energy (reference) PCA profile (2.02- 
2.82 kcV), yielded the time lags as measured with ISGRI 
and HEXTE above 20 keV. Note, that a time shift of 
-t-0.5 X 125/iS has been applied to the PCA time stamps, 
because the PCA times refer to the start of the time pixel 
instead of the mid in the case of the E_125us_64M_0_ls 
event mode. Both the PCA/HEXTE and ISGRI time lags 
as a function of energy are shown in Fig. El covering 
the 2-100 keV energy band. We see that the ISGRI and 
PCA/HEXTE time lag measurements in the overlapping 
20-35 keV energy band are fully consistent, and that the 
time lags decrease for increasing energies, attaining values 
consistent with zero for energies in excess of 35 keV con- 
sistent with a trend already set in near ^ 6 keV. The full 
consistency of the PCA/HEXTE and ISGRI time lags in 
the overlapping 20-35 keV band indicates that the time 
alignment between RXTE and INTEGRAL is better than 
~ 25/xs. 
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Fig. 5. Pulse-profile collage of IGR J00291+5934 using data from RXTE/PCA (2-58 keV; panels a~e), RXTE/SEXTE 
(20-151 keV; panels g-j), INTEGRAL/mU-X (5-10 keV; panel f) and INTEGRAL/lSGRl (20-150 keV; panels k- 
n). Two cycles are shown for clarity. All profiles have high off-set values due to the nature of the instruments. The 
error bars represent 1 sigma statistical errors. All profiles reach their maximum near phase ~ 0.95. Notice the highly 
sinusoidal shape of the profiles for energies up to 100 keV. 



3.3.3. Spectrum of the pulsed emission and pulsed 
fraction 

Important diagnostic parameters for constraining the 
parameter space in theo retical modelling (see e.g. 
IViironen fc PoutanenlfioO^ are the pulsed spectrum and 
the pulsed fraction (here defined as: pulsed fiux/total flux) 
as a function of energy. While the pulsed fluxes can be 
derived for all the main instruments considered in this 
study, the total flux of IGR J0029 1-^5934 can only be 
derived reliably for the high-energy instruments with ar- 
cminute imaging capabilities. In particular, the nearby 
source V709 Gas contaminates the total flux estimate of 
IGR J00291-h5934 using the non-imaging RXTE instru- 



ments (see section r^. 2 .1(1 . Galloway et al. (2005) did use 
the PGA to derive the fractional rms amplitude of the 
pulsations as a function of energy (in their Fig. 2, top 
panel); therefore, their values contain a contribution from 
V709 Gas and cannot be used to estimates the fractional 
rms amplitude of IGR J00291+5934. 

We derived the pulsed fluxes for the PGA, HEXTE and 
ISGRI instruments (the JEM-X statistics for the pulsed 
component is too low, cf. Fig. O panel f ) and the total 
fluxes only for the imaging INTEGRAL instruments'^. To 
derive the pulsed fluxes we determine the number of pulsed 



^ The RXTE and INTEGRAL measurements were nearly 
contemporaneous. 
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Fig. 6. Time lags as a function of energy in the 2- 
100 keV energy range combining RXTE/PCA (2-35 
keV; 2004 December 7-10; open circles), RXTE/REXTF. 
(20-60 keV; 2004 December 7-10; filled diamonds) and 
INTEGRAL/ISGRl (20-100 keV; December 3.6-9.8; filled 
squares) measurements. 



excess counts in a given energy band. A truncated Fourier- 
series fit using 3 harmonics to the pulsed phase distribu- 
tions as a function of energy yielded this quantity. In case 
of the PCA the pulsed excess counts (2-30 keV) have been 
converted to flux values (ph/cm'^skeV) in a forward en- 
ergy spectral folding method assumin g an underlying ab- 
sorbed (iVn = 2.8 x lO^^ cm'^; see lNowak et al.ll2004l) 
PL model with an energy dependent index taking into ac- 
count the different exposures and energy response of the 
5 different PCUs. Figure shows these PCA pulsed flux 
measurements in an EF^ representation. 

For HEXTE an equivalent method - dividing the 
pulsed excess count rates in a certain energy band by its 
effective sensitive area assuming a PL model - was used 
to obtain the pulsed fluxes in the 15.6-267.5 keV energy 
range taking into account the different cluster and 1 
dead time corrected ON exposures and energy responses 
(see Fig.[7I). 

Finally, for ISGRI the pulsed excess counts (20-300 
keV), accumulated during an exposure of 336.94 ks, have 
been converted to flux values using (a) the pulsed excess 
counts of the Crab pulsar as a reference from an equiv- 
alent 5x5 on-axis dither observation (Rev-102; 54.33 
ks), determined in the same energy bands as chosen for 
IGR J00291+5934, and (b) a spectral model for t he pulsed 
emission of the Crab pulsar (see e.g. iKuiper et al...2001i|) 
over the 15-10 000 keV energy range based on HEXTE and 
CGiJO/COMPTEL measurements. These ISGRI pulsed 
flux measurements are also shown in Fig. 13 

From the PCA/HEXTE and ISGRI pulsed flux mea- 
surements shown in Fig. 13 it is clear that this emis- 
sion component is hardening towards higher energies. 



E [kcV] 

Fig. 7. The total and pulsed spectrum of 
IGR J00291-I-5934 in an EFe representation from 
3 to 300 keV combining flux measurements from 
i?Xr£;/PCA/HEXTE and INTEGRAL/lSGRl for the 
pulsed part and INTEGRAL/ISGRl and JEM-X for the 
total part. Notice the hardening of the pulsed spectrum 
towards higher energies. 

Superposing now also the total IGR J00291-I-5934, spec- 
trum from JEM-X and ISGRI measurements (see section 
13.2. l|l in Fig. [7\ which softens weakly towards higher en- 
ergies, it is clear that the pulsed fraction increases as a 
function of energy. The energy dependence of the pulsed 
fraction is shown in Fig. |H1 It gradually increases from 
- 6% at 6 keV to - 12 - 20% near 100 keV. It is the first 
time that such behaviour has been measured from any of 
the known accretion-powered ms pulsa r s. 

We note here that Gallowa v et al.l l)2005l) report de- 
crease of the fractional rms amplitude of oscillations 
with increasing energy above 10 keV. As remerked be- 
for, this can be explained by the contribution in the PCA 
data from V709 Gas, which has a harder spectrum than 
IGR J0029H-5934 in this energy range. 

4. Discussion 

4.1. Companion star and the distance 



Taking the INTEGRAL observed total fluence of 1.37 x 
10""^ erg cm~^, and a recurrence time T of 3 years 
ljR,emillardll2004l) . we find a long term time-averaged flux 
of 1.45x10"-^^ erg cm~^ s~^ (smal ler than the 1.8xl0~^ ^ 
erg cm^-^ s^^ flux reported bv iGallowav et all l2005l) 
Following the curr ent model ofVerbunt fc van den Heuvel 
(1995) and Bildst en fc Chakrabartv. (^200l|) . we therefore 
find, for an inclination i = 90° and minimum companion 
mass Mc = 0.039 M0 (for a 1.4 Mq neutron star), a lower 
limit on the source distance of 4.7 kpc. At a most proba- 
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Fig. 8. Tlie pulsed fraction (=pulsed flux/total flux) 
of IGR J0029H-5934 using pulsed/total flux measure- 
ments from PCA/HEXTE and JEM-X/ISGRI (open cir- 
cles, PCA pulsed flux relative to JEM-X total flux; open 
squares, HEXTE pulsed flux relative to ISGRI total flux; 
filled squares, ISGRI pulsed flux relative to ISGRI total 
flux). The pulsed fraction gradually increases from ~ 6% 
at ~ 6 keV to ~ 12 - 20% at ~ 100 keV. 

ble inclination angle i = 60°, we obtain a source distance 
of 5.4 kpc. In both cases the companion is a hot brown 
dwarf. The ~1 kpc difference betw een our calculated min- 
imum source distance and that of iGallowav et al"! ^|2005^ 
is attributable to the difference in observed total flux. At 
a distance of - 5 kpc, IGR J0029H-5934 lies about 270 
pc off the galactic plane, placing it within the thick disk, 
but outside the thin disk. 



4.2. Origin of X-ray emission 

The spectral parameters obtained from our ther- 
mal Comptonization fits to the broad-band data for 
IGR J00291-H5934, kT^ = 49 kcV and tt = 1.12, are sim- 
ilar to those of other MSPs: kT^ = 60 keV, tt = 0.88 in 
SAX J18 08.4-3658, fcT, = 33 keV, tt = 1.7 in XTE J1751- 
350 fsee iGierlihski fc Poutanenll2005(l . We note that the 
product Tt x kT^ is amazingly close in the three sources. 
The spectral shape and the product tt x kT^ are very sta- 
ble during th e outburst as is obser ved in other sour ces too 
l|Gilfanov et al. 1998.: _Gicrlihski fc Poutanenll2005() . 

The constancy of the spectral slopes during the out- 
bursts and their extreme similarity in different MSPs can 
be used as an argument that the emission region geom- 
etry does not depend on the accretion rate. If the en- 
ergy dissipation takes place in a hot shock, while the 
cooling of the electrons (that emit X/7-rays via ther- 
mal Comptonization) is determined by the reprocessing 
of the hard X-ray radiation at the neutro n star surface 
(so called two-phase model, see lHaardt fc MaraschL1993t 



Stern et alJll995t IPoutimCTlll998t iTitarchuk et alJll998l: 



Malzac et al.ll200lll . the spectral slope is determined by 



the energy balance in the hot phase and is a function of 
the geometry. At constant geometry (e.g. slab), the tem- 
perature depends on the optical depth, but tt x kT^ is 
approximately constant. 

4.3. Pulsed fraction and pulse profile 

We showed that the pulsed fraction gradually increases 
with energy from ~ 6% at 6 keV to ~ 12 — 20% near 100 
keV. For the first time such behaviour has been measured 
for any of the accretion-powered MSPs. 

A spot emitting as a black body (with isotropic spe- 
cific intensity) at a slowly rotating star pr oduces nearly 
sinus oidal variations in the observed fiux I Beloborodovl 
and the peak is reached when the spot is closest to 
the observer. When the spot radiation is anisotropic, the 
harmonics appear in the pulse profile, wi th their strength 
depending on the degree o f anisotropy ijPoutanenI l2004t 
IViironen & Poutanenll20o3) . In case of rapid rotation, as 
in MSPs, Doppler boosting affects the variability pat- 
tern. When the spot moves towards the observer, the 
emission increases, while for a spot moving away, the 
flux drops. The Doppler factor reaches the maximum a 
quarter of the period before the peak of the projected 
area shifting the emission peak towards earlier phase. 
The observed flux due to the Dop pler effect varies as the 
Doppler factor in power (3 -1- F) l|Poutanen fc Gierlihskil 
I2OO3I: IViironen fc PoutcmCTll2004(l . where the photon in- 
dex F could be a function of energy. If the Doppler factor 
varies around 1 with 2% amplitude, we get 10% variability. 

The Comptonized spectrum can be approximated as 

(1) 



FEO,E-(^"-'')exp{-[E/E,f) 



where Ec ~ kT^ is the energy of the cutoff and parameter 
/3 ~ 2 describes its sharpness. The local photon index is 
then 



TiE) EE 1 - 



dlnFe 
dln£; 



^To + f^iE/E,) 







(2) 



At low energies, F « Fq, and rms (or pulsed fraction) is 
a very weak function of energy. Close to the cutoff, the 
spectral index rapidly increases and the pulsed fraction 
should grows with energy, as observed (see Fig. ^ . The 
Comptonization models predict softening of the total spec- 
trum with simultaneous hardening of the pulsed spectrum 
at higher energies. 

The observed pulse profi les are almost purely sinu- 
soidal (see Fig. |5l and also ICallowav et al.l l2005t) . The 
total variability amplitude is proportional to the prod- 
uct sin i sin 9 fPoutanc n 200 4) , where 9 is the angu- 
lar distance of the spot from rotational axis. The rela- 
tive amplitude of the second harmonics to the fundamen- 
tal is also proportional to the same product. The ~ 6% 
obs erved amplitude constrains the angles i and 9 (see 
e.g. iGierlinski fc Poutanenll2005l) . We estimate that for a 
I.4M0 neutron star with 12 km radius, 9 varies between 
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3° and 15° when the indmation i varies from 90° to 
18°. 

4.4. Time lags 

The observed spectrum of IGR 100291+5934 consists of 
a black body from the neutron star surface and a compo- 
nent produced by Comptonization of these seed photons 
in the hot electron region, presumably a shock which can 
be represented as a plane-parallel slab. The angular distri- 
butions of the black body and Comptonized photons emit- 
ted by the slab are significantly different. The difference in 
the emission patterns causes the two components to show 
a different variability pattern as a function of the pulsar 
phase, with the hard Compto nized component leading the 
soft black body component IjPoutanen fc GierlinskilEoO^t 
IViironeu PoutauenlEioi. 

This scenario is consistent with observations of 
SAX J1808.4-3658 llCui et all Il998l: iGierlihski et all 
I 20021) . XTE ,11751-350 llGierlihski fc PoutanenI l20n5r 
where time lags increase rapidly with energy until 7-10 
keV, where the black body contribution becomes negligi- 
ble and then saturate above 10 keV. In IGR J00291-h5934, 
however, time lags first increase until 7 keV and the n 
decrease slightly (see Fig. and iGallowav et al.l l2005l) . 
Judging from the PCA data, it seems that the lags sat- 
urate above 15 keV. The decrease of the time lags above 
7 keV does not have an obvious explanation. There are 
also (not very significant) indications from the HEXTE 
and ISGRI data that the lags even reach zero value at 
~ 50 keV. If confirmed, this kind of behaviour would be a 
serious challenge to any model. 

5. Summary 

We analyzed the spectral and the timing behaviour from 
the entire INTEGRAL observation of IGR J00291-h5934, 
which was spatially well distinguished from the neighbor- 
ing source V709 Gas. The time averaged broad-band spec- 
trum from 5 to 300 keV obtained with JEM-X/ISGRI is 
well described by a thermal Comptonization model with 
seed photons from the neutron star surface scattered in 
a shock-heated accretion column above the hot spot re- 
gion. Additionally, studying the spectral evolution, we 
found that the electron temperature and the related opti- 
cal depth arc almost invariant as the flux decreases during 
the outburst. 

The INTEGRAL/lSGRl data, supported by simulta- 
neously obtained RXTE/PCA and HEXTE data, also al- 
lowed us to study the timing behaviour in the millisec- 
ond range for energies up to 150 keV. We detected the 
1.67 ms period consistently with all three instruments and 
found that the pulse profile is close to a sine function 
at all stu died energies as observed for the other MSPs 
(see e.g., iGierlinski et al.l I 20021 : ICierlinski &: PoutanenI 
l2005|l . This is the first time that pulsed emission has 
been detected from 2 up to ~ 150 keV for any of the 
currently known accretion-powered millisecond pulsars. 



Moreover, it demonstrat es the excellent tim ing capabili- 
ties of INTEGRAL fsee lKuioer et al.ll2003t) down to the 
0.1 ms level. We have discovered that the pulsed fraction 
significantly increases with energy which is naturally ex- 
plained by the effect of Doppler boosting on the exponen- 
tially cutoff Comptonization spectrum. We have measured 
a spin-up rate of -1-8. 4x 10~^^ Hz s~^. Accurate measure- 
ments of the frequency derivative are necessary to con- 
strain models of N S spin evolution in X-ray binaries (e.g., 
lT.a,mb fc Yull2004^ . 

From our total fluence measurements assuming that 
the mass accretion rate is driven by gravitational radi- 
ation, we find the minimum source distance of 4.7 kpc. 
This locates the source, like other accreting X-ray MSPs, 
within the thick disk of the Galaxy. 
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